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Abstract 

In this work, we study the emission of tensor-type gravitational degrees of freedom from 
a higher-dimensional, simply rotating black hole in the bulk. The decoupled radial part 
of the corresponding field equation is first solved analytically in the limit of low-energy 
emitted particles and low-angular momentum of the black hole in order to derive the ab- 
sorption probability. Both the angular and radial equations are then solved numerically, 
and the comparison of the analytical and numerical results shows a very good agreement 
in the low and intermediate energy regimes. By using our exact, numerical results we 
compute the energy and angular-momentum emission rates and their dependence on the 
spacetime parameters such as the number of additional spacelike dimensions and the 
angular momentum of the black hole. Particular care is given to the convergence of our 
results in terms of the number of modes taken into account in the calculation, and the 
multiplicity of graviton tensor modes that correspond to the same angular-momentum 
numbers. 



1 Introduction 



It has been more than a decade since the introduction of the new theories postulating the 
existence of additional spacelike dimensions in nature; namely, the large extra dimensions 
[1] and warped extra dimensions scenarios [2] have led to intense research activity of 
the theoretical as well as the phenomenological consequences of that existence. The 
introduction of a new, lower than the four-dimensional, fundamental scale for gravity has 
created the expectation that the elusive quantum theory of gravity might manifest itself 
soon during high-energy particle collisions, with energies higher than the new gravity 
scale that may be as low as a few TeVs, at ground-based accelerators. The products 
of these collisions will inevitably be manifestations of a strong gravity theory, including 
possibly the creation of higher-dimensional miniature black holes P]. 

These equally elusive objects, although very small, will be created and decay in 
a controlled environment and in front of our detectors with the emission of Hawking 
radiation [1] being their most distinctive feature. As a result, there has been a consid- 
erable amount of interest in the study of the radiation emission spectra from a higher- 
dimensional black hole in the literature in recent years (for some reviews, see O [6]). 
Before reaching the Plank phase, a miniature black hole passes consequently through 
balding [7], spin-down, and Schwarzs child phases [SIE])- The starting point was the class 
of spherically symmetric black holes that were considered to describe the longest and 
more important Schwarzschild phase in the life of these objects. The corresponding emis- 
sion of Hawking radiation was exhaustively studied, both analytically and numerically 
0|9l[IOl[IIl[T2l[T3l[Tll[l5l[T6l[n], leading to a number of interesting results. 

However, the most generic type of a black hole produced by the collision of two 
particles with a nonzero impact parameter is a rotating black hole. As a result, the 
interest was eventually turned to the study of the axially symmetric spin-down phase 
[l8l[T9l|20linil22l|23l|2l|25l|26l|271^ that was initially considered to be significantly 
shorter than and preceding the Schwarzschild phase. However, a recent Monte Carlo 
simulation [29] (see also [30]) that has included the effect of rotation of the black hole, 
has found that the "spin-down" phase is not as short-lived as it was thought and that a 
separate Schwarzschild phase with no angular momentum might not exist at all. 

The study of the spin-down phase is important for an additional reason: the ques- 
tion of the energy balance [31] between the "bulk" and "brane" channel during the 
Hawking radiation has not been answered yet. Studies of the Schwarzschild phase 
[9l [ini [131 El 113 [ini [U] that included the emission of both scalar fields and gravi- 
tons have revealed that the brane channel is in most cases the dominant one, although 
at certain circumstances the bulk channel can be equally important at specific particle 
channels. Similar studies have also been performed in more recent years [321 1321 [Ml [2S] 
for the emission of scalar fields during the spin-down phase, with the dominance of the 
brane channel still persisting. 

In order to give a final answer to the energy balance question, we also need to in- 
clude in our calculations the emission of gravitons during the spin-down phase. Until 
recently, the field equations of gravitational perturbations in a higher-dimensional, ax- 
ially symmetric black-hole background were not known. Even today, we have at our 
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disposal the field equations of specific gravitational modes in certain classes of axially 
symmetric gravitational backgrounds. The perturbation equations for tensorlike gravi- 
tational modes in the case of a higher- dimensional rotating black hole with D > 7 and 
equal angular-momentum components was derived in [36]. In a subsequent work, one of 
the authors of this work derived the corresponding equations for tensor-type gravitons 
for higher- dimensional black holes with one angular-momentum component and D > 7 
again [23 • Then the stability and quasinormal modes of the considered tensor-type grav- 
itational perturbations were investigated in [38]. Recently, in [2S] perturbation equations 
were derived for particular scalar, vector and tensor-type gravitational modes for a five- 
dimensional rotating black hole with two equal angular-momentum components. 

In this work, we study the emission of Hawking radiation in the bulk in the form 
of tensor-type gravitational modes by a higher-dimensional black hole with one angular- 
momentum component. The geometrical background will be the one considered in P?] . 
and we will therefore demand the existence of at least three additional spacelike dimen- 
sions. Our analysis starts in Sec. 2 where we present in more detail the geometrical setup 
and theoretical framework. In Sec. 3, we solve analytically the radial part of the graviton 
field equation in the limit of low-energy emitted particles and low angular momentum 
of the black hole, and we derive an analytic expression for the absorption probability. 
We then move, in Sec. 4, to solving the set of angular and radial equations via numer- 
ical methods. Our numerical techniques are presented in Sec. 4.1 with our results for 
the absorption probability described in Sec. 4.2. In the same section, we also perform a 
comparison of the analytical and numerical results to check the validity of the analytic, 
approximate method. The exact form of the energy and angular-momentum emission 
rates in the bulk from the simply rotating black hole in the form of tensor-type gravitons 
are finally computed and presented in Sec. 4.3. We finish with the presentation of our 
conclusions in Sec. 5. 



2 Theoretical framework 

As mentioned in the Sec. 1, the perturbation equations for gravitons in a higher- 
dimensional, rotating black-hole background have been derived in a limited number of 
cases. Here, we will focus on the case considered in where the higher- dimensional 
gravitational background can be described by a line element of the form 

ds^ = Gmn dz^^ dz^ = gab dx'' dx^ + S\x) dVtl, (1) 

where {a, b} = (0, 1, 2, 3) and c/fi^ stands for the line element of an n-dimensional unit 
sphere 5". The above line-element is a special (4 -(- n)-dimensional case of a more general 
class of gravitational backgrounds where the spacetime can be written as the warped 
product of an m-dimensional spacetime N' and an n-dimensional space JC of constant 
curvature [IHl E]. In these type of backgrounds, gravitational perturbations can be 
classified into tensor, vector and scalar types according to their transformation properties 
as tensors on the const ant- curvature spacetime /C. It was this property that allowed for 
the derivation of perturbation equations for all types of gravitational modes in the case 
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where 



of a maximally symmetric higher-dimensional black-hole background [12] where m = 2 
and /C = 5". 

The Myers-Perry solution [43J that describes a D-dimensional black hole with = 
[{D — l)/2] independent angular-momentum parameters does not, in general, belong to 
the aforementioned class of line elements. However, when the black hole rotates only 
in a single two-plane along the four-dimensional spacetime, its line element takes the 
well-known form 

+ (^r' + a' + ^-^^^ym' 9 d^' + r' cos' 9 dnl (2) 

A = r' + a'-J^, j: = r' + a' cos'9, (3) 

and where we have set ai = a and a2 = = . . . = = 0. The above assumption 
is justified by the demand that the black hole is created by the collision of particles 
propagating on our four- dimensional brane: the colliding particles have a nonzero impact 
parameter only along the usual 3-space which results in a single angular-momentum 
component once the black hole is created. Then, the black hole's mass Mbh and angular 
momentum J are related to the parameters a and /i as follows |43] : 

MsH=^-^^^j^^l^, J=^M,^a, (4) 
IottGd n + 2 

with Gd being the (4 + n)-dimensional Newton constant, and An+2 the area of an (n + 2)- 
dimensional unit sphere given by A„+2 = 27r*^"+'^^/^/r[(ri + 3)/2]. 

The line element ([2]) is a special case of the class of backgrounds considered in 
Jini SI], and more specifically of the class described by Eq. ([T]) with S{x) = rcos9. 
As first stated in [27] and later demonstrated in more detail in [H], the tensor-type 
gravitational perturbations for the line element ([2]) - which exist only for n > 3, or 
D > 7 [37] - can be expanded in terms of a basis of transverse and traceless harmonic 
tensors T^^'"'* on the unit sphere as follows: 



SG,, = 2S\x) , (5) 

where {i,j} refer to the y coordinates along the sphere S'", and T-^'"'' satisfy the eigen- 
value equation 

[A + £(£ + n-l)-2]Tg'") = 0. (6) 

In the above, A is the Laplace- Beltrami operator on 5", and £ = 2, 3, 4, ... an integer 
number that labels the corresponding eigenvalues. Finally, a is a label to distinguish 
harmonic tensors with the same eigenvalue. 
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Under the expansion the component of Einstein's equation in vacuum leads 
to the following second-order hyperbolic equation for the amplitude h!^'°'\x) 

- BHt - g^' da{r cos 6) d,HT + t ^ Ht = 0, (7) 

r cos u cos"^ U 

where □ is the d'Alembertian operator for the metric gab{x), and where, for simplicity, 
we have omitted the labels {i,a}. Under the further factorization 

Ht{x) = e"^' e™^ R{r) Q{e) , (8) 

the above partial differential equation reduces to a set of radial and angular equation, 
namely, 

r" V A r 



d,(r-AdrR)+( — - ' -A,,^)i? = 0, (9) 



■ / „^ de (sin e cos" 9 dgQ) + (u^a^ cos^ 9 - ^ - ^ , + E,em) Q = 0. 
sm^cos"^ V sm^^^ cos"' 6* / 

(10) 

In the above, we have used the definitions 

K = [r + a )uo — am , rs.jim = -E-jfm + o w — 2amu , (11) 

with Ejim being the separation constant of the two equations, and j a new quantum 
number that labels the eigenvalues of the angular function Q{9). 

In order to obtain the complete solution for the wave function of the tensor-type 
gravitational perturbations of the background ([2]) one needs to solve the above set of 
second-order ordinary differential equations (l9l)-(fT0l) for R and Q. The same is, in princi- 
ple, necessary for the computation of the energy emission rate for Hawking radiation in 
the form of tensor-type gravitational degrees of freedom: the radial equation will yield the 
expression for the absorption probability (or graybody factor) for the particular type of 
particles, with the angular equation providing the value of the separation constant Ejim 
that appears in the former equation. The above task can be performed either analytically, 
in the low-energy and low-angular- momentum limit, or numerically with no restriction 
on these two parameters. In the next two sections, we will follow both approaches to 
fulfill this task. 



3 Analytic solution 

As was noted before [23 SI] in the context of more general analyses, when the spacetime 
background has the form of Eq. ([1]), the tensor- type gravitational perturbations are 
found to satisfy the same field equations that a massless scalar field obeys in the same 
background. In the present case, the same result also holds as the set of equations (l9l)- (fT0l) 
are identical to the ones that follow from the scalar field equation 

-i= dM (v^ G*^^9^<l>) = , (12) 
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if the following expansion of $ in terms of the hyperspherical harmonics Y^^'") (y) on 5*" 
is used 



-iu}t 



(13) 



Given the different nature of the scalar and gravitational degrees of freedom, the two 
sets of equations differ only in the allowed values of the angular-momentum number i: 
whereas, in the scalar case, it satisfies the constraint £ > [46J, this changes to £ > 2 
in the case of gravitons ^7\. The decoupled set of equations for a massless scalar field 
propagating in the higher-dimensional background ([2]) first appeared in and were 
further used in [331 EH] for the study of the energy emission rates for Hawking radiation 
emitted by the simply rotating Myers-Perry black hole in the form of scalar fields in the 
bulk. 

The radial equation IQ was analytically solved for scalar fields propagating in the 
bulk in [33| in the low-energy and low- angular-momentum approximation. As the equa- 
tion for tensor-type gravitons is identical, apart from the allowed range of values for 
i, the analytic solution in this case follows along the same lines. For this reason, here 
we give only a brief account of the analysis and the results obtained for the absorption 
probability in the analytic approach, results that are compared with the exact numerical 
ones in the next section. 

The analytic approach amounts to finding first the asymptotic solutions near the 
horizon of the black hole (r ~ r^), and far away from it (r ^ r^) and matching the two 
in an intermediate zone, to create an analytical solution for R{r) over the whole radial 
regime. The black hole's horizon radius Vh follows from the equation A{rh) = 0, and is 
given by the relation r^"*"^ = yu/(l + a^), where a* = a/vh- 

In the near-horizon regime (r ^ Vh), Eq. (Q can be rewritten in the form [33] 



/(l-/)^ + (l-/^./)f + 



in terms of the new radial variable [27] /(r) = A(r)/(r^ + a^) and the quantities 



(n + 1) + (n - l)at 



(1 + a^)a;* - a^m . 



R = 0, 
(14) 

(15) 



In the above, we have also defined = urt and = 1 — Aal/A^. Equation f|T^ can be 
brought to the form of a hypergeometric differential equation with the general solution 



RNH{f) = A^ r{l-f)^F{a,b,c; f) 

+ A+r-(l-/)'^F(a-c+l,6-c+l,2-c;/) 



(16) 



where A± are integration constants, and the indices (a, b, c) are defined as a = a -|- /3 + 
D^, — 1, b = a + /3, and c = 1 -|- 2a. In addition, the parameters a and (3 are given by 
a± = ±iK^/A^ and 



/3 



K2-[i(i + n-l)al + A,f,m]{l + al 



{D, - 2)2 - 4 



Al 



(17) 
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Close to the horizon, the general solution f lTB]) can be written as the sum of an incoming 
and an outgoing plane wave after employing a convenient transformation of the radial 
variable, namely, y = r/i(l + al) ln(/)/A*. If we impose the boundary condition that no 
outgoing modes exist near the black hole's horizon, we can set either = or = 0, 
depending on the choice for the sign of a. The two choices are found to be equivalent, 
thus we choose a = «_ and A^ = 0. Then, the near-horizon solution acquires the form 

RNHif) = A. r (1 - ff F{a, b, c; f) . (18) 



On the other hand, in the far-field regime (r ^ r^j), Eq. (|9]) can easily be brought 
[33] into the form of a Bessel differential equation if we make the substitution R{r) = 

r~^~^^ R{r) and employ a new radial variable z = ur, 
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If we further define for convenience the quantity = y Ejim + o?^"^ + (^^) ) the general 
solution in the far-field regime may be written as 

RpFir) = ^Ju K) + ^Y, {ur) , (20) 

where Jy and Yy are the Bessel functions of the first and second kind, respectively. 

Before the two asymptotic solutions (JT8l) and ( l20l) can be matched, they both need 
to be expanded for intermediate values of the radial variable. To this end, the hypergeo- 
metric function appearing in f lTSj) needs also to be shifted so that its argument changes 
from / to 1 — / by using a well-known relation |33] . Then in the limit r ^ rh or, 
equivalently / — )■ 1, the near- horizon solution takes the "stretched" form 

RMr) ~ r -("+1) ^ + A,r ^^^'^ , (21) 

with Al and A2 defined as 

,/3 r(c)r(c — a — b) 



Al = A-[(l + aDrn' 



T{c-a)T{c-b) ' 



A - A rn I n^w"+ir^^+''*~'^ r(c)r(a + &-c) 

A2 - A_[(l + ajr, J r(a)r(6) " ^^^^ 

Similarly, the far-field solution (|20|) in the limit of r — )■ takes a similar polynomial form 

r 2 1 (z/ + Ij TT r 2 l^-^- j 

For the two stretched solutions to perfectly match, the power coefficients of r need to be 
the same. It can be easily shown that this is indeed the case in the limit of a* < 1 and 
<1. Then, by ignoring terms of order (w^, a^, a^u^) or higher in the expressions of /3 
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and z/, we find that (ra + 1) /3 ~ — j, (n + 1) {/3 + — ~ — (j + n + 1) and z/ ~ j + 
By identifying the coefficients of the same powers of r, we finally obtain the constraint 

Bi ^ (2/a;rfe)^^^+"+^ u T\u) T{a + f3 + D,- 1) T{a + (3) T{2 - 2/3 - D,) 

T-> 2)+n + l 5 l^^J 

B2 vr (1 + a2)^iT^ r(2/3 + - 2) r(2 + a - /3 - D,) r(l + a - P) 

that guarantees the existence of a smooth, analytic solution for the radial part of the 
tensor-type graviton wave function for all r, valid for small and w*. 

Let us, at this point, clarify the expression for the eigenvalue Eji^ that appears 
both in /3 and v. This quantity does not exist in closed form, and can be found either 
numerically or in terms of a power series in the limit of small auj. We reserve the use of 
the first method for the next section - in the context of the present analysis, valid in the 
low-energy and low-angular-momentum, we may use instead the analytic power series 
expansion [HI H9] 



E,i^ = Y.fk {auf . (25) 

fc=0 

For the accuracy of our analysis, we keep terms up to 4th order - the exact expression^ 
of the coefficients fk can be found in |18]. It is only in the expansion of the power 
coefficients of r in the matching process that all terms beyond the first one are ignored; 
in this case, Ejim — fo = jU + n + 1), where j > i + \m\ and ^ jg^ Z+}. 

A quantity that determines, to a great extent, the Hawking radiation emission rate 
of the black hole is the absorption probability \Ajern\'^ - or graybody factor, since it is the 
reason for the deviation of the black-hole spectrum from a pure blackbody one. We may 
derive it, by expanding the far- field solution (|T9l) for r — )• oo, in which case we obtain 



RMr) 



n + 2 



r 2 v27ra; 



(Si +i52)e-H-^-f '^-fJ + (5i -iE2)eH'^^-t'^-f) . (26) 



The absorption probability is then easily determined via the amplitudes of the outgoing 
and incoming spherical waves, namely. 



Bi - iB. 



Bi + iB2 



2 



2i {B* - B) 



BB* + i{B* - B) + l 



(27) 



where B = B1/B2 is given by Eq. 1^^. The above result can be used to evaluate the 
absorption probability for the emission of tensor-type gravitons in the bulk, from a simply 
rotating black hole, in the low-energy and low-angular- momentum limit. 



4 Numerical analysis 

In this section, we use numerical analysis in order to solve both the angular and radial 
equations for any value of the energy of the emitted particles and angular momentum 

^For consistency, we should point out that in [A'Sl the indices are interchanged compared to the 
ones in this work, and the total sign of the /2 coefBcient should be reversed due to a typographical error. 
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of the black hole. We start by presenting the main aspects of our numerical techniques, 
and then we turn to the derivation of exact numerical results for the graybody factor 
and energy emission spectrum for tensor-type gravitational modes in the bulk. 



4.1 Numerical techniques 

The angular equation ffTOj) . in terms of the new variable x = cos(26'), can be written in 
the following form 

2(1 - x^) Q"{x) + [n-l-{n + 3)x] Q'{x) (28) 

+ — ^ + + 7 7— Qi^) = . 

\ 2 4 X — I 1 + X J 

The above differential equation has been solved in the literature before in different con- 
texts and forms: for instance, its four- dimensional version (that follows for i = n = 0), 
in the presence of a positive cosmological constant, was solved in [50]; in |18], the above 
higher-dimensional version was solved for scalar fields, i.e., for £ = 0, 1, 2, . . .; finally, for 
tensor-type gravitons living in a higher- dimensional space but in the presence of a neg- 
ative cosmological constant, the corresponding equation was numerically solved in [44] . 
Thus, the numerical analysis demanded for solving Eq. (|28l) for tensor-type gravitons 
living in a higher- dimensional asymptotically flat spacetime is a simplified case of the 
one presented in [33] - we refer the interested reader to that work for more details. 

The differential equation (!28|) has three regular singular points, at x = ±1 and 
X = oo. The angular function Q{x) can be alternatively written as 

g(x) = (l-x)H/2(i + a;p^(;,), (29) 

where x = 2^; — 1. Under further expansion of the rescaled function y{z) in terms of an 
infinite series of Jacobi polynomials, supplemented by regularity conditions aX z = and 
z = 1, Eq. fl28|) takes the form of an algebraic equation - a three-term recurrence relation 



Cfc+i Ok + Ck(3k + Ck-i 7fc = (30) 

for the coefficients Ck appearing in the expansion of y{z). In the above relation, the coef- 
ficients ak, (3k, and 7^ are constants depending on the fundamental parameters (a;*, a*, n) 
of the theory, the quantum numbers {i,m), the angular eigenvalue Ejem, and the new 
index A; = 0, 1, 2, . . . that labels the power of the expansion. 

The solution for the angular function Q{6) is of limited physical importance for the 
calculation of the energy emission spectra for the simply rotating black hole. On the other 
hand, the computation of the eigenvalue Eji^, that also appears in the radial equation 
dH]) as a separation constant, is of paramount importance. The value of the separation 
constant can be obtained by applying the infinite continued fractions method pT]. The 
continued fraction equation [44J follows from the three-term recurrence relation fl30|) and 
involves ratios of successive terms of the coefficients ak, Pk, and 7^. This equation can 
be numerically solved in any desired accuracy for the value of Ej£m, for given values of 
00^,, a^,, and n. 
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In the case of vanishing angular momentum of the black hole, the value of the sepa- 
ration constant can be found in a closed form by the requirement that the corresponding 
power series expansion of y{z) with a finite number of terms converges [51]. In that case, 
we obtain [HI Hi] 

(31) 



Ejirn = {2k + i+ \m\){2k + i + \m\ + n + 1) = j{j + n + 1) 



where, in the last part of the above equation, we have set j = 2k + i + \m\. Under this 
alternative definition, the eigenvalue is labeled by a new quantum number with values 
j = 2, 3, 4, . . . and, at its lowest order, it coincides with the one for a (n + 2) sphere, in 
agreement with the discussion below Eq. (123]) . 

If the rotation parameter a of the black hole is nonvanishing, the eigenvalues Ejim are 
in principle noninteger and complex. In that case, we can find the value of Ejim{co*, a^,), 
for any value of and a^:, by using the following procedure. 

1. We start from the nonrotating black hole and find the exact value of Ejim, for the 
corresponding j, according to Eq. ( 13T1) . 

2. We increase the rotation parameter by a very small amount and search for the 
closest to the previously found solution for Ejim- 

3. We repeat the previous step until any required value of a* is reached and all cor- 
responding values of Ejim are found. 

By following the aforementioned process, we are able to compute the values of the angular 
separation constant Ejim, for any and a*, and thus to proceed to the numerical 
integration of the radial equation (Q. 

Equation (Q can in turn be rewritten in an alternative form under the redefinition 
of the radial function R{r] 
coordinate defined through the relation dr^ 
reads 



r "/^ (r^ + a^) P(r) and the employment of the tortoise 



a?)dr//S.. The new equation then 



rf2p(r, 
dr'i 



+ 



am 



u — 



+ 



A 



(r^ -|- a'^y 



U{r) 



Pin 



(32) 



where 



f/(r) 



+ n - l)a^ 



+ A 



nin 



2) 



4^2 



+ 



3a2 



+ 



(n + 1)/. _ 



1 



n 

2^ 



+ 



+ a^ 



(33) 



In this form it is straightforward to derive the asymptotic solutions at the horizon and 
spatial infinity. First, at the horizon, if we set r — )■ r/^ and A — t- 0, we easily obtain 



Pin) ~ e'^"* + A2 e 
where ^41^2 are integration constants, and 

bJ = bj — m^h = liJ 



am 



a^ + r 



2 ' 



(34) 



(35) 
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with fi/i the rotation velocity of the black hole. Since no outgoing wave is allowed to 
classically exist outside the horizon of the black hole, the physically relevant solution of 
Eq. ( 132|) at the horizon is 



-i(u -mVLh) -T 



rl + 



P{r) = A2 e-^^'-* ~ (r - r,) ^'(^'^) {Z^ + 0{r - r,)) , (36) 

where is a rescaled integration constant. As expected, the near-horizon solution f[T5]) . 
derived in Sec. 3, reduces to the same expression if we take the limit / — )■ 0, expand A 
in powers of (r — r^) and redefine the integration constant A_. 

For the purpose of our numerical analysis, we introduce close to the horizon the new 
function 

,H . (1 _ ^)-''''*"'""''"Vm. (37) 

Since P(r) satisfies the asymptotic condition ( l36l) . z{r) is regular at the event horizon. 
We may also fix the value of any undetermined integration constant, by setting 

z{rn) = 1. (38) 

If we then expand z{r) near the event horizon as 

z{r) = 1 + z'ivh) (r -rh) + 0{r- r^f , (39) 

and substitute into Eq. fl32|) . we find the value of z'{rh) which, together with Eq. (138|1 . 
are the boundary conditions for our Eq. (132!) at the horizon. 

Next, at spatial infinity (r — )■ 00), the two linearly independent solutions of Eq. fl32|) 

are 

P,(r)~e-^'^^ Po(r)~e*'^^ (40) 

which describe the ingoing and outgoing wave, respectively. The functions Pi{r) and 
Po(r) can be found analytically as series expansions for large r up to any order. 

The numerical integration of Eq. ( l32l) then proceeds as follows: with the eigenvalue 
Ejirn already numerically known for all values of and a*, we start from the horizon, 
with the values of z{rh) and z'{rh) as boundary conditions, and move outwards by using 
the ND Solve built-in function in Mathematica@ for < r < rf, where Vf ^ r^. After 
the function P(r) is known numerically, we find a fit of this function by considering the 
superposition of the two solutions f l40p in some region near r^: 

Pir) = Z,P,,{r) + Z,P,{r). (41) 

The fitting procedure allows us to find the coefficients Zi and Zo. In order to check the 
precision of the coefficients we increase the internal precision of NDSolve, the value of 
Tf, and the number of terms in the series expansion for Pj(r) and Po(r), making sure that 
the values of Zi and Zq do not change within the desired precision. The same shooting 
procedure, though for different boundary conditions, has been used recently in for 
analysis of stability of higher- dimensional black holes. 
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Figure 1: (a) Comparison between our analytical (solid lines) and numerical (data points) 
results for the graybody factor for the modes (j = £ = 2, m = 0) and (j = 5, £ = 2, m = 
1), for a = 0.5 and D = 7. (b) Absorption probabilities for the sets of modes with 
i = 2, from left to right, (j = 2,m = 0) (red ine), (j = 3, m = —1,1) (green ines), 
(j = 4,m = —2,0,2) (blue lines), and (j = 5,m = —3,-1,1,3) (magenta lines), for 
a = 0.5 and D = 7. 



Once this process is completed, the quantity Zo/Zi gives the ratio of the amplitudes 
of the outgoing and ingoing modes at a large distance from the black hole, and the 
absorption probability follows easily through the relation 

\Ajim\'^ = 1 — \T^jim\'^ = 1 — \Zo/Zi\'^. 



4.2 Absorption probability 

By following the two approaches described in Sees. 3 and 4.1, we have derived analytical 
approximate results as well as exact numerical ones for the absorption probability for 
gravitational tensor modes that propagate in the background of a higher- dimensional 
simply rotating black hole. The two sets of results ought to agree in the low-energy and 
low- angular-momentum limit, but we expect them to deviate once we move outside these 
regimes. In order to check the extent of the agreement of the two sets of results as well as 
its dependence on the particular mode studied, in Fig. [T]^a) we depict these two sets for 
two indicative modes with (j = 2, £ = 2, m = 0) and (j = 5, £ = 2, m = 1): the analytical 
results are given by the solid lines whereas the numerical results are presented as data 
points - both sets of results correspond to the case with D = 7 (or n = 3) and a = 0.5 
(in units of r/j). As expected, the agreement between the two sets is indeed very good in 
the low-energy and even intermediate-energy regime, but inevitably it breaks down as 
we move towards the high-energy one. The agreement is better for the lowest modes and 
it worsens for higher modes for which the graybody factor raises to a significant value 
and approaches unity at an increasingly higher value of the energy parameter uvh- 

In Fig. Wljo), we examine the aforementioned behavior of the graybody factors for 
different tensor modes by using exact numerical results. We classify the modes primarily 
by the angular-momentum number j which can be considered as the total angular- 
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Figure 2: Absorption probabilities for the mode (j = 2,£ = 2,m = 0) as (a) a function 
of n = 3, 4, 5, 6, 7, for a = 0.5, and (b) a function of a = 0, 0.5, 1, 1.5, 2, 2.5, for n = 3. 



momentum number of the mode, with £ denoting the angular-momentum along the 
compact space and m the one in the plane of rotation of the black hole. As in the 
case of scalar fields [33l[35], a set of modes corresponds to each value of j: the constraints 
j > i + \m\ and ^ |o, Z"*"} |18] dictate that for each value of j, i can take values 

in the range [2, j] while, for given j and i, m can take j — i + 1 values in total. In Fig. 
[D^b), we display the set of modes corresponding to the values j = 2, 3, 4, 5 - in order to 
keep the plot tidy, we fix £ = 2 and present the graybody factors for the modes with 
the j — i + 1 allowed values of m in each case. We may clearly see that as either j or m 
increases, the corresponding graybody curve shifts to the right and to higher-energies - 
a similar behavior would have been observed if we also varied i. 

Next, we investigate the dependence of the gravitational tensorial graybody factors 
on the spacetime parameters of the theory, namely the number of additional spacelike 
dimensions n and the angular- momentum parameter of the black hole a. In Fig. |2]^a), we 
display the absorption probabilities for the indicative mode [j = 2, i = 2,m = 0) as n 
changes from 3 to 7, while keeping the angular-momentum parameter fixed at a = 0.5. 
The graybody factors for the gravitational modes in the bulk clearly decrease as the 
number of transverse-to-the-brane spacelike dimensions increases. For the same mode, in 
Fig. [2t^b), we present the dependence of the graybody factors as a changes from to 2.5, 
while keeping the dimensionality of spacetime fixed at = 7. In this case, the graybody 
factors for tensorlike gravitons are clearly enhanced as the angular-momentum of the 
black hole increases. This behavior is in total agreement with the one observed for bulk 
scalar fields [331 ES] propagating in the same background. 

4.3 Energy and angular- momentum emission rates 

Having determined the exact value of the absorption probability, we can now proceed to 
compute the differential emission rates of energy and angular momentum from a higher- 
dimensional simply rotating black hole in the bulk in the form of tensor-type gravitons. 
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These are given by the following expressions, 



dujdt 

j,e,m 



^j:^^;j^l^sASn\A,em\^ (42) 



dudt 2n 



^E^s7^^^T(^")iA.™r, (43) 



where u is defined in Eq. f p5|) and the temperature Th of the black hole is 

{n + 1) + {n - l)al 



47r(l + a2)r,, 



(44) 



The quantity Ngrp is the multiplicity of the second-rank symmetric, traceless {T^a = 
0) and divergence-free {DbT^^ = 0) tensor harmonics Tab that satisfy Eq. ([6]). Equiv- 
alently, it is the multiplicity of tensor modes on S*" that, under the aforementioned 
constraints, are described by the same angular-momentum number i. This number was 
calculated by Rubin and Ordonez in [52] and found to be 

, _ {n + l)(n - 2){n + £)(/ - l)(n + 2£ - l)(n + £ - 3)! 

^s^^^ ' - 2(£+l)!(n-l)! ' ^^^^ 

for the £th eigenvalue. The above formula was derived by expanding the tensor harmonics 
Tab in terms of the harmonic functions Y"^^^ and utilizing the representation theory 
of SO{n + 1). In the Appendix, we give an alternative proof of the formula (l45ll . in 
which we start from a similar expansion for the tensor harmonics Tab and determine the 
multiplicity of modes by finding the number of independent solutions of the equations 
that the expansion coefficients satisfy - these equations express the properties of Tab on 
S"', such as its symmetry, its tracelessness, and transversality. This method is completely 
elementary and can also be extended into antisymmetric tensor harmonics and higher- 
rank tensor harmonics. 

In order to compute the differential rates fl42p - (H5]) . we need to sum the contribution 
of all tensor modes labeled by the different values of the (j, i, m) angular quantum 
numbers. In practice, the sums need to be truncated at an appropriate high value of 
each number in such a way that the derived values of the two rates are as close as 
possible to the real ones. To this end, we adopt the following procedure: we first fix 
one of the angular numbers and sum over the other two within this range - in this 
way we find the contribution of each value of the fixed parameter to the total sum. If 
the contribution of the highest considered value of an angular number is not small, we 
increase the particular value range. We repeat the described procedure for all angular 
numbers until the contribution to the energy and angular momentum emission rates of 
the highest considered multipole number becomes negligibly small. 

As an indicative example, in Figs. |3t^a,b) and IH^a,b) we display the contributions of 
the lowest m and I tensor modes to the energy and angular-momentum emission rates, 
respectively, for D = 7 and a = 0.5. In all cases, we may observe the increasingly smaller 
contribution of the higher modes to the specific rate, and thus the convergence of the 
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Figure 3: The energy emission of tensor-type gravitons for D = 7, = 0.5 (black line) 
together with the contributions of different quantum numbers (color lines) are shown. In 
the left-hand figure: m = —2 (cyan line), m = — 1 (blue line), m = (green line), m = 1 
(red line), m = 2 (purple line), m = 3 (magenta line), m = 4 (orange line), m = 5 (brown 
line), m = Q (gray line). In the right-hand figure: i = 2 (green line), i = 3 (orange line), 
i = 4 (red line), i = 5 (magenta line). 

The largest contribution correspond to m = 0. Peaks of positive m contributions (m = 
1, 2, 3, 4, 5, 6) lay to the right from the peak of m = contribution. 



corresponding sum. This is due to the fact that, according to Fig. [T]^b), the higher modes 
become important at a larger value of the energy parameter ur^, and this in practice 
takes place after the peak of the emission curves - determined by the temperature of 
the black hole - has been reached. As a result the higher modes contribute mostly to 
the "tail" of the emission curves. This is more clearly shown in Fig. [5] where the energy 
emission rate is presented, for D = 7 and a = 1, in terms of the contribution of the j 
modes: as the highest considered value of j increases, from j = 5 to j = 8, then to j = 10, 
j = 12, and finally to j = 15, the emission curve becomes wider and the slope of the tail 
decreases, whereas the low-energy behavior and the peak of the curve remain unchanged. 
As mentioned above, in all the cases studied in this work, care was taken so that the 
change in the emission curves would be negligibly small when a cutoff was imposed on 
the highest values of all angular numbers. In general, as either n or a* increases, the 
number of modes that need to be summed increases, too - in order to obtain as accurate 
as possible emission spectra, we have summed up to j = 22, i.e., m = 20, in the cases 
considered. 

Next, we turn to the dependence of the energy and angular-momentum emission 
rates of the black hole on the spacetime parameters, namely, n and a^,. In Figs. ISIa,b), 
we illustrate the dependence of the energy spectrum on the number of additional spacelike 
dimensions and the angular momentum of the black hole, respectively. As in the case 
of scalar field^ [33| |35] . the energy emission rate has a very strong dependence on n 
with an enhancement of almost two orders of magnitude as n changes from n = 3 

^As a check of our numerical analysis, we have successfully reproduced the exact results for the 
energy emission rate of scalar fields in the bulk from a higher-dimensional simply rotating black hole 
[55] that were derived with an independent code. 
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Figure 4: As above, for the angular-momentum emission of tensor-type gravitons (black 
line). In the left-hand figure: m = —2 (cyan line), m = —1 (blue line) m = (green 
line), m = 1 (red line), m = 2 (purple line), m = 3 (magenta line), m = 4 (orange line), 
m = 5 (brown line), m = 6 (gray line). In the right-hand figure (from top to bottom): 
£ = 2 (green line), i = 3 (orange line), i = 4 (red line), i = 5 (magenta line). 
The largest contribution corresponds to m = 1. Peaks of other positive m contributions 
(m = 2, 3, 4, 5, 6) lay to the right from the peak of m = 1 contribution. The contributions 
of negative values of m are negative. 
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to 77. = 7. This enhancement is present in all energy regimes with the emission curve 
becoming significantly taller and wider as n increases. The dependence of the energy 
spectrum on the angular momentum of the black hole is also nontrivial, although of a 
smaller magnitude: for the case n = 3 depicted in Fig. Mjo), the increase of the angular- 
momentum parameter from a* = to a* = 1.2 results into the emission of less energy 
per unit time in the low and intermediate regime and an enhancement in the emission of 
high-energy modes. As n gets larger, this dependence becomes milder, a feature which 
is again in accordance with the behavior of the bulk scalar fields emitted by the same 
black hole spacetime. 

In Figs. [71^a,b), we depict the dependence of the angular-momentum emission rate 
on the same spacelike parameters. As the number of extra dimensions increases, we 
observe again a significant enhancement in the rate of loss of angular momentum by 
the black hole. This enhancement reaches more than an order of magnitude and results 
in the emission of a higher number of modes in all energy regimes. Contrary to what 
happens in the energy spectrum, the increase in the rotation velocity of the black hole also 
increases the angular-momentum emission rate from the black hole. The enhancement is 
significant, although of a smaller magnitude than the one in terms of n, leads to the loss 
of angular momentum via the increased emission of modes in the whole energy spectrum, 
and manifests itself independently of the dimensionality of spacetime. 

Finally, in Fig. |8]one can see the total energy emission and angular- momentum emis- 
sion for the tensor- type gravitons for D = 9 and a fixed a* = 1.2 as well as contributions 
of different quantum numbers m, calculated by the accurate shooting method. In Table 
I the total emission power by the scalar field for a^, = 1.0 (taken from J35j) is given in 
comparison with the total emission power of tensor-type gravitons. There we can see that 
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Figure 5: The energy emission of tensor-type gravitons for D = 7, a* = 1. The different 
(upper) curves, from left to right, correspond to the highest value of j considered in the 
sum: j = 6 (red line), j = 8 (green line), j = 10 (blue line), j = 12 (magenta line), and 
J = 15 (black line). 



n 


Scalar field 


Tensor-type gravitons 




3 


0.1646 


0.0013 


0.8% 


4 


0.3808 


0.0222 


5.8% 


5 


0.7709 


0.1853 


24% 



Table 1: Total emission power (mass loss rate, in units of 1/r^) by scalar field (a* = 1.0 
taken from [35j) and by tensor-type gravitons (a* = 1.2). 

although at small number of spacetime dimensions n the contribution of gravitons into 
the total radiation is tiny, it quickly increases with n and becomes dominant for large n. 



5 Conclusions 



Because of the complexity of the analysis demanded for the derivation of the perturbed 
gravitational equations in the background of a higher-dimensional non-spherically- 
symmetric black hole, the emission of Hawking radiation in the form of gravitational 
modes from such a spacetime has been up to now an uncharted territory. In this work, we 
have used the results of a previous analysis performed by one of the authors [37] according 
to which the derivation of the field equations for tensor-type gravitational perturbations 
is indeed possible under the assumption that the spacetime manifold is the warped prod- 
uct of two submanifolds with its line element having the form of Eq. (1). This class of 
spacetimes includes not only the previously studied higher- dimensional spherically sym- 
metric black-hole backgrounds but also the case of a (4 -|- n)-dimensional rotating black 
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Figure 6: The energy emission rate for tensor-type gravitons in the bulk for (a) a^, = 1 
and n = 3, 4, 5, 6, 7 (from bottom to top), and for (b) tt, = 3 and a^, = (blue line, top), 
= 0.5 (green line), a* = 1.2 (red line, bottom). 



d^J/dt doj d~}/dtdu, 




Figure 7: The angular-momentum emission rate for tensor-type gravitons in the bulk for 
(a) a=K = 1 and n = 3,4,5,6,7 (from bottom to top), and for (b) n = 3 and a* = 0.5 
(green line, bottom), a^: = 1.2 (red line, top). 



hole with a single angular-momentum component along the (3-|-l)-dimensional brane. In 
addition, it was further shown that these equations, upon the use of tensor harmonics 
as a basis, can lead to a set of decoupled ordinary differential equations with respect to 
the spacetime coordinates. 

The derived equations for the tensor-type gravitons propagating in the bulk are found 
to be identical in form with the ones satisfied by bulk scalar fields. We were thus able 
to analytically study the problem of the computation of the absorption probability by 
using techniques employed previously for the emission of scalar fields by the same type 
of black hole, under proper modifications to allow for the different values of the angular- 
momentum quantum numbers that characterize the graviton modes. This study led to an 
analytical expression for the absorption probability for tensor-type bulk gravitons valid 
in the limit of low-energy emitted modes and low-angular-momentum of the black hole. 

In order to derive the complete emission spectra, for arbitrary values of the energy 
of the emitted mode and angular momentum of the black hole, we also performed an 
exact numerical analysis to solve both the angular and radial part of the graviton's 
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Figure 8: The energy emission (left) and angular- momentum emission (right) of tensor- 
type gravitons for D = 9, = 1.2 (black line) together with the contributions of different 
quantum numbers m (color lines) are shown. 



field equation. In the process, the value of the angular eigenvalue, that appears in and 
connects the two equations and which does not exist in closed form, was also computed. 
Having all the above exact results at our disposal, we were thus able to find the value 
of the absorption probability, or graybody factor, for tensor-type graviton modes in the 
specific background. 

The exact form of the absorption probability in terms of the energy parameter ur^ 
was studied in detail as well as its dependence on the particular graviton mode considered. 
A comparison between the approximate analytical and exact numerical results for its 
value was performed, and it was found that, for the lowest graviton modes, the agreement 
of the two sets of results is remarkably good and extends up to the high-energy regime; 
as higher modes are considered, the analytical result deviates from the exact one at 
an increasingly smaller value of the energy. The dependence of the graybody factor 
on the spacetime parameters, namely the number of additional spacelike dimensions and 
angular momentum of the black hole, was also investigated. According to our results, the 
absorption probability for tensor-type gravitons decreases with the number of transverse- 
to-the-brane dimensions but increases as the black hole rotates faster - this behavior is 
similar to the one found for bulk scalar fields in previous analyses [331 ES] • 

We next moved to the computation of the emission spectra, namely the energy 
and angular-momentum ones. The value of the corresponding differential emission rates 
strongly depends on the number of modes - characterized by the set of (j, i, m) angular- 
momentum numbers - that are considered in the sum. For this reason, we performed 
a careful study of the convergence of our results before imposing a cutoff on the three 
quantum numbers. In all cases studied, we made sure that the effect of all the higher 
modes left out of the sum was always negligible. In addition, a technical calculation was 
performed for the derivation of the multiplicity of tensor modes characterized by the 
same set of angular-momentum numbers - this number is distinctly different from the 
one for bulk scalar fields and affects the value of the differential emission rates. 



Combining the above, the energy and angular-momentum emission spectra were 
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finally computed. Both spectra exhibit a very strong dependence on the number of ad- 
ditional spacelike dimensions with the increase in the rate of emission of either energy 
or angular momentum reaching even 2 or 1 orders of magnitude, respectively. The de- 
pendence on the angular momentum of the black hole is more particular: while the 
angular-momentum emission is clearly enhanced, the differential energy emission rate 
displays either an enhancement of the high-energy modes and suppression of the other 
frequencies, for low values of n, or a rather mild dependence of the spectrum on a^,, for 
high values of n. 

Apart from its obvious theoretical interest, the calculation of the emission spectra of 
a higher-dimensional, simply rotating black hole in the form of gravitons in the bulk has a 
very important phenomenological interest in the exciting case of the creation of miniature 
black holes at ground-based accelerators. Previous studies [33l [35] have revealed that 
the bulk emission of the other species of particles allowed to propagate in the whole 
spacetime, namely, the scalar fields, is subdominant compared to the emission that takes 
place in the form of brane- localized scalar fields. When this is combined with the fact that 
the total number of fermionic and gauge bosonic degrees of freedom of Standard Model 
are also restricted, and thus emitted, on the brane, the brane emission channel becomes 
even more dominant. Addressing the question of energy balance between the brane and 
bulk channel for the last species, i.e., the gravitons, is of paramount importance for the 
estimate of the percentage of the total energy of the black hole which is lost in the 
bulk, and thus of the chances for the potential detection of the produced black hole via 
the emitted Hawking radiation. The analysis and results presented in this work on the 
emission of tensor-type graviton modes in the bulk are the first necessary step towards 
this direction that will hopefully be soon complemented by a similar calculation for the 
vector and scalar gravitational modes. 
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Note added. While this manuscript was being written, a complementary work [S3] 
appeared that also studies the emission of tensor-type gravitons in the bulk by a simply 
rotating black hole. 

A Multiplicity Formula for Tensor Harmonics on 

The multiplicity on S*" of the second-rank symmetric tensor harmonics Tab, satisfying 
Eq. dS]) along with the conditions T^a = and DbT^"^ = 0, was calculated by Rubin 
and Ordonez |52j . They started from the expression for an eigen harmonic tensor 

T,, := TAsD^n^.D^n,, = ^ fefl J F(L)(fi) , (A.l) 

£(m) 
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in terms of the harmonic functions Y"^^^ with the eigenvalue l{l + n — 1), represented 
as polynomials of the Cartesian coordinates f2* for the Euclidean space E^^^ in which 
5" is embedded as O • O = 1. Here, Da is the covariant derivative with respect to 
the standard metric 7ab = ■ of S*". This implies that Ti^i^ transforms as the 
product of two irreducible representations of SO(n + 1), (2, 0, ■ ■ ■ ) x (£, 0, ■ ■ ■ ). Then, 
they showed that among the irreducible representations obtained by decomposing this 
product representation, only the (^, 2, 0, ■ ■ ■ ) has the same eigenvalue as Tj^j^ for the 
standard send-order Casimir operator and its multiplicity is one. Therefore, iVgrp(S'"') 
coincides with the dimension of the irreducible representation (£, 2, 0, ■ • • ), leading to the 
formula ( 145|) . 

In this Appendix, we give an alternative proof of the formula f H5|) . which as we 
mentioned in Sec. 4 is more powerful as it can be extended into antisymmetric tensor 
harmonics and higher-rank tensor harmonics. The starting point is the following expres- 
sion for Tab similar to flA.ip : 

Tab = a,,...j^.,,,,^V^ ■ ■ ■ 9.^'Da^''Db^'\ (A.2) 

where aj-^...ji,-i-^i^ is a set of constants satisfying the following conditions: 

^ji---je',iii2 ^iji---je)\hi2y (A. 3a) 

^ji---je-2j'' ;hi2 ~ 0' (A. 3b) 

Hh-ji-Ji+i)' = 0- (A-3c) 

Utilizing the formula DaDb^-' = —■Jab^'', it is easy to see that Tab defined by (lA.2p 
satisfies the eigenvalue equation (E]), and the conditions on are equivalent to the 
conditions T^a = and D^T^^ = 0. 

What we would like to do is to count the independent number of solutions to the linear 
algebraic conditions ( 1A.3p . The basic idea is to reduce the problem to lower dimensional 
problems by decomposing the range of the indices i,j,--- = 1, ■ ■ ■ , n -|- 1 to jl = + 1 
and ■ ■ ■ = 1, ■ ■ ■ , n. 

We start from the condition (IA.3cp . This condition reads 

P 

for z = tl and 



P 



+ l-p 

for i = i'. Hence, the condition = can be rewritten as 



(A.5) 



a^-m = 0' ««-si':tt« = 0' «s-tt:r = 0' «ttoi-:i;i;+i) = O' ^{ji-fi-r.+di' = (A.6a) 

«S-ttjl-:i^S = (2<P<^+1), (A.6b) 
p(^p — 1) 

""hM-r.-M = ^^^i_p)^f^^2-pf ^---^^''^--<-^''-^ (2 < P < ^). (A.6c) 
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Next, we consider the condition flA.3b|) . For ii = i[ and ^2 = ^^^^ condition can 
be divided into the following set of conditions: 



(A.7) 



Hence, taking into account the expressions ( ]A.6ap . we find that all components of aj^...-ij^i 
can be expressed in terms of ajj...j^:j'^j^ and satisfying the constraints 



as 



T\\tzP] k 



«tt-tt;i'tt = 

p{p — 1 



2 Ja 



(0<p<£-2), 

(2 <J9<^+1), 



(£+1 - p)(£ + 2-p) 



(A.8) 

(A.9a) 
(A.9b) 

(A.9c) 

(A.9d) 

(A.9e) 



Similarly, for ii = ^ and 12 = i', the condition ( lA.3bl) can be decomposed into the 
sequence of conditions 

(^t-U-r,-'i + ""i-U'rr-r.n = o (o < p < ^ - 2), (A.io) 

which can be written with the help of ( ]A.9aP as 

+ + 1 - p)^mi,.,..,^ = (0 < p < ^ - 1). (A.ii) 

Because this sequence of equations is consistent with contraction, if they hold for p = 
i — l,i — 2, they hold for the other values of p as well. Further, the equation for p = i — 1 
is a contraction of the first of ( 1A.8I) . Hence, the only independent condition is 



0. 



(A.12) 



Finally, for ii = ^2 = tl^ the condition (]A.3bp can be decomposed to 

H-ij[-%m + H-U'3[-r,.n = o (o < p < £ - 2). (A.is) 

Inserting the relation (lA.9aP into these, we get 

(p + l)(p + 2)(£ + 1 - + 2 - P)am%...ik'k'f,...-,'^_,r.) = (0 < p < £ - 2). 

Thus, the number of degrees of freedom of a* is identical to the number of linearly 
independent solutions to the constraints f lA.Sp . f lA.12p and flA.14p for ctjj^.. o'o' and 
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We first calculate the degrees of freedom of ajj...j^:i'^j^. Let us introduce the following 
symbols representing tensors: 





K 


k' 




-2m + l-Je-2m+2) ^ 


= 0, 


■■■Am]), 


(A. 14a) 






k' 
' ' k' 


1 k' 

-1 {jl'"if_2m- 


^l-h-2m+2l'^rn 


= 1, 


■■,[(^ + l)/2]), 


(A. 14b) 




k[ 

%{■ 


k' 
' ' k' 


k' 




I,--- 


, [^/2] + 1) 


(A.14c) 




K 

«4' 


k' 
' ' k' 


"-^ '^m-l'^m ., 
-2 ■'l 


Je-2m + 2 m-l'^m 


m = 


2,---,[^/2] + l). 


(A.14d) 



In terms of these symbols, the conditions on ajj...j^.j'^j^ are given by A^^^ in flA.SP and 



5^™) : (4m-3)(£-2m + 2)(£-2m+l)X'" + 2m(2m-l){2(£-2m + 2)F'" + Z™} 
(m = 2,3,--- ,[£/2]). (A 

First, we show by induction that under the condition A^^\ B^'^^ holds for any m = 
2,3,--- , [^/2] if B^"^^ holds. Let us assume that holds. Then, from the above 

contraction formulas, the following condition holds: 

(E^"^))' : (4m-7)(£-2m + 2)(£-2m + l)X" + 4(2m2-m-4)(£-2m + 2)y'" 
+ (4m2 - 2m - 8)Z"' + 8(2m - 3)(m - 1)W"' = 0. (A.16) 

Further, from the condition A^^\ for each level m, we obtain 

A^"^^ : {i-2m + l)X'^ + 2mY"' = 0, (A.lTa) 
4""^ : {i - 2m + 2)Y"' + Z"" + 2{m ~ 1)W"' = 0. (A.17b) 

It is easy to find that these conditions are related by 

^(^M)' _ 4(2m - 3)v4S™^ - 2(2m - 3){i -2m + 2)A'^'^^ = (m - 2)B^"'\ (A. 18) 

This implies that i?*^™) holds. Hence, we only have to impose the condition B'^'^\ 

Now, we calculate the dimension of the linear set A*^^^ With the notations 
j'li - ■ ■ ,i" = 1, - ■ ■ ,n — 1, y4*^°) is equivalent with 

fln-ninn = 0, Pan-niJ'-i^Yi-p:"" = - P + l)«n-nOY-:j;Yi-p)" (P = 1' • • • ' (A. 20) 

«n---n:m" = 0, PO-n-nj'-/ ■■■j'^_^_j^_^:ni" = (i — P + '^)0'n-n(j'{-:j';^^^^)i" (P = 1, . . . , ^). (A. 21) 

Hence, the number x„ of linearly independent equations of A^'^^ obeys the following 
recurrence relation: 

n+e-iCe+1 + nn+i-iCi, Xi — 1 (A. 22) 

This can be solved to yield x„ = Un+iCi+i. 
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Next, we count the linearly independent equations of B^'^\ Let us rewrite B^"^^ as 

12annf^...f^'_^.,r,n = " " " > 2), (A.23a) 
^Qannnj'^-j'/_.^:nn = " " " > 3), (A.23b) 

{p - l){5p - Q)a^...nj-...f;_^:nn = " " " > 4, p = 4, ■ ■ ■ , £) , (A.23c) 

where the right-hand sides of these are linear combinations of terms of the form 
an...nj'-^...j'^'_^:nn < p) and a...,j'i/'. Thcse equations are obviously linearly independent, 
and the number of them is n+e-3Ce~2- Thus, because we can show that A^*^^ and B^'^^ are 
independent, the number of degrees of freedom of ctjj - j^ii'ji^ is given by 

= ^ 2 — ~n+e-iCe — nn^^Ci^i — n+i-?,Ce,-2- (A. 24) 

The degrees of freedom of Ottir-if i:«i«2 counted in the same way. As in the 

previous case, let us introduce the symbols X^", and Wp that are obtained 

from the definitions ( ]A.14al) by inserting jj at the top position of the indices as XJ" : 

• ■ -tZ Ui-jL2m-i^jL2mj^-2^+i) ("^ = 0, " " " , [{i-l)/2]). Thcu, iu tcrms of these symbols, 
the conditions on a-jjj - j^ jii'jii, are given by 

: {£ - l)4k'(j[...r,_,.f,_,y + ^' ji---r,_,f,_,.k'i' = 0, (A.25a) 

Blf"^ : {Am-l){i-2m + l){i-2m)Xp + 2m{2m + l) {2{i-2m + l)Yf' + Z^} =0 
(m = 0,l,-- - ,[(£+l)/2]). (A.25b) 

First, we show by induction that under the condition A\^\ if B^^^ is satisfied, B^^^ is 

also satisfied for any m = 1, 3, ■ ■ ■ , [£/2]. Let us suppose that B^"^~^^ holds. Then, from 
the contraction rules above, we obtain 

: (4m-5)(£-2m)(£-2m + l)Xj" + 4(2m2 + m-4)(£-2m + l)yj"^ 
+ (4m2 + 2m - 8)Z™ + 8(2m - l)(m - l)Wp. (A.26) 

Further, applying the contraction rules to the condition A^^^ leads to the following two 
equations at each level m: 

A'f^ : {£-2m)X^ + {2m + l)Yf = 0, (A.27a) 
A'^^^ : {i-2m + l)Yp + Z^ + {2m-l)Wp = 0. (A.27b) 

We can show that these equations are related by 

(2m + - 8(m - - 8(m -!)(£- 2m + 1)^5"^^ = (2m - 3)5j'"\ (A.28) 

This implies that B^"^^ holds. Hence, we only have to impose the condition B^^\ 

Next, we count the linearly independent equations of the constraints. First, the con- 
dition Ap^^ can be written 

(P + ^Wfl-f;_^:m' = ■■■ , (p = 2, ■ ■ ■ , £), (A.29) 
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where the right-hand side of this equation is a hnear combination of terms of the form 
^n---nj'-^---j'/_^:nn (<? < p) and a...;jii'i. Hence, the number of hnearly independent equations 

is nn+e-3Ci_2- Next, the number of independent equations of B^^^^ is obviously n+eCi+i. 
Hence, because we can show that A^^^ and B^^^ are independent, the number of the total 
degrees of freedom of CLij[-j'^_^:i[i'^ is 

Ni^ = ^ ^ — ^-n+e-2C£-i — n+iCe+i — nn+e~3C£-2- (A. 30) 

Finally, adding the numbers of the degrees of freedom of and «jtjj---j;_-^:i;ii,, we 

obtain 

Ngrp(S"') = ^ ^ — ^ {n+e^iCi + n+e-2Ce-i) — {n + 1) {n+fCe-i + n+i-3Ci-2) ■, (A. 31) 
which after some algebraic manipulation reduces to Eq. f l45|) . 
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